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The observation of quantum oscillations [l|, in 
underdoped cuprates has generated intense de- 
bate about the nature of the field-induced re- 
sistive state and its implications for the 'normal 
state' of high Tc superconductors. Quantum oscil- 
lations suggest an underlying Fermi liquid state 
at high magnetic fields H and low temperatures, 
in contrast with the high-temperature, zero-field 
pseudogap state seen in spectroscopy. Recent 
heat capacity measurements [3] show quantum os- 
cillations together with a large and singular field- 
dependent suppression of the electronic density of 
states (DOS), which suggests a resistive state that 
is affected by the d-wave superconducting gap. 
We present a theoretical analysis of the electronic 
excitations in a vortex-liquid state, with short 
range pairing correlations in space and time, that 
is able to reconcile these seemingly contradictory 
observations. We show that phase fluctuations 
lead to large suppression of the DOS that goes 
like a/T? at low fields, in addition to quantum os- 
cillations with a period determined by a Fermi 
surface reconstructed by a competing order pa- 
rameter. 

The 'normal state' of the high Tc superconducting 
cuprates remains an enigma. In the underdoped regime, 
close to the Mott insulator, a large pseudogap (~ 50 
meV) is observed in the electronic excitations by a va- 
riety of spectroscopic and thermodynamic probes for 
T > Tc y, [sj in the absence of a magnetic field H. It 
thus came as a great surprise when quantum oscillations 
were observed 0, Q in the low-T regime, once supercon- 
ductivity is destroyed hy H > iJiir, the irreversibility 
field. Such oscillations, periodic in 1/H, are most easily 
understood in terms of a Fermi liquid state. This raises 
the question: what is the role of the strong correlation 
physics that leads to the high-T, zero-field pseudogap in 
the low-T, high-field quantum oscillations? After all a 50 
T field cannot destroy a 50 meV pseudogap. 

A very important clue comes from recent electronic 
specific heat measurements 0] which see quantum oscil- 
lations riding on top of a strongly suppressed DOS with a 
singularity at low fields. The unusual 7?-dependence 
points toward a d-wave gap node-like structure @ per- 
sisting even into the resistive state. It is not a priori clear 
how to reconcile the Fermi surface probed by quantum 
oscillations with the d-wave gap. 



Further, the large suppression of the specific heat 7 ~ 5 
mJ / molK^ &i H = 50T (compared with the normal state 
value 7 ~ 18 mJ/molK2 at = 0) implies that the one 
has not recovered the 'normal state' at the experimentally 
accessible fields. This is consistent with observations of 
nonlinear diamagnetism Q and large Nernst efl[ect Q 
seen in the normal state and suggestive of substantial 
phase fluctuations above Tc- There is also considerable 
evidence that the superconducting transition in under 
doped cuprates is associated with phase-disordering [13- 
121 rather than a gap collapse. 

In this paper we make a simple model of the electronic 
excitations in a vortex liquid and show how it helps rec- 
oncile, within a single framework, the apparently con- 
tradictory observations described above. In the vortex 
liquid state the local d-wave pairing amplitude is non- 
zero, but superconducting correlations are short-ranged 
in both space and time. Our analysis generalizes earlier 
studies of the mixed state of s-wave superconductors in 
refs. [H-Iil in two ways - d-wave pairing and dynamical 
phase fluctuations - both of which are very important 
for quantum oscillations in cuprates. We show that the 
effect of phase fluctuations on the electronic self-energy 
leads to quantum oscillations riding on top of a strongly 
suppressed DOS that goes like ^fH at small H . 

However, we need more than just phase fluctuations to 
understand the quantum oscillations. Their observed fre- 
quency in underdoped cuprates (unlike that in overdoped 
samples [l^) is known to be too small to be consistent 
with a Luttinger Fermi surface (FS). It corresponds to 
an electron-like FS with area of only about 2% of the 
Brillouin zone (BZ) [H, [l]. It is now accepted that these 
observations necessarily imply a FS that has been re- 
constructed by a (possibly field-induced) density wave 
order 
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Thus, to get a complete description of 
the underdoped cuprate experiments, we incorporate both 
phase fiuctuations and a competing order parameter at 
the end of the paper. 

Phase Fluctuations: Wc characterize the vor- 
tex liquid state with a simple ansatz for the 
gauge- invariant correlation function D^i^(r,t) = 
(vI/^(r,t)M/:(O,0)exp(z|f/„^A.dl)). The field ^^ir,t) 
describes singlet pairs on the bond (r, r -|- afi.) with a 
being the Cu-Cu lattice spacing of the Cu02 square 
lattice with sites r; /i = ztx, zty and A is the vector 
potential for the magnetic field H = Hz. The retarded 
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FIG. 1: Quantum Oscillations in Density of States: DOS at the chemical potential Nh{0) normalized by the zero-field 
normal state A'^o(O). The magnetic field H, measured in units of /x, is huc/jJ. with ujc = eH/m*c. We use A ~ 0.3, and an 
impurity broadening 70 = 0.01 in units of = 1. a: Nh{0) vs. 1/H for the static and dynamic cases (F 7^ 0) analyzed in the 
Landau levels (LL) basis. The DOS for LL's with A = is also shown. Note the (l/H) periodicity of the oscillations and the 
difference in the damping between the static and dynamic cases. The inset shows the self-energy approximation (see text), b: 
iZ-dependent suppression of DOS computed in the LL basis. We show that Nh{0) oc for small H. c: DOS obtained from 
the k-space and semiclassical quantization schemes for F — 0.1vf/£. 



correlation function 

(r, t) = s^^A' cxp {-ry2f) cxp (-rt)e(t) (1) 

is assumed to be short ranged in space and time; Q{t) is 
the Heaviside step function. A is the local d-wave pair- 
ing amplitude that persists above i?irr, as expected for a 
superconductor where the resistive transition is governed 
by phase fluctuations. The d-wave nature is described 
by s^iy = 1 for /i = dzv and —1 otherwise. The spatial 
decay in ([T]) is on the magnetic length scale £ = y/hc/eH, 
set by the average inter- vortex separation in the extreme 
type-II limit. We work in a regime where the cyclotron 
radius i?c ^ ^ ^ Co ^ ^f^' where is the vortex core 
radius and kp^ the interparticle spacing. 

The temporal decay in ([1]) is governed by an energy 
scale HT that characterizes vortex motion. On general 
grounds we expect < T < vp/i- The upper limit 
arises from ballistic motion of vortices with the Fermi 
velocity vp. For simplicity, we write T = avp/i with 
< a < 1 in the dynamic case where vortices are mo- 
bile. The separable form ([!]) simplifies the algebra, but a 
more elaborate non-separable correlator with T = Dq^, 
where D is the vortex diffusion coefficient, is not ex- 
pected to change our conclusions qualitatively (2^, [2^ . 
We will also find it useful to compare our results for dy- 
namic phase fluctuations (F 7^ 0) with the the static case 



' /2£ ) with time-independent 



Dfj,t,{r,t) = s^^A^exr 
phase fluctuations [l3Hl5j 

We use the simplest approximation for self-energy 
E(k, oj) (inset of FiglT]a) to find the effect of phase fluctu- 
ations on electronic excitations in the vortex-liquid state. 
Our approach generalizes the static, s-wave analysis of 
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refs. IJ, A similar self-energy has also been used for 



the pscudogap phase of cuprates [22 - 24 1 , but no calcula- 
tions have been presented for quantum oscillations. 



The central quantity of interest to understand quan- 
tum oscillations is the single-particle DOS at the chemical 
potential iV(0) at T = as a function of the external field 
H. We use the self energy S to compute the electronic 
Green's function G = I/[G(7^ — E], where Go is the free 
Green's function. The imaginary part of G then gives us 
the DOS N{ui). We note that there is no anomalous part 
of the Green's function, since (\E'p(r,t)) = in the ab- 
sence of long-range phase coherence. We will focus first 
on the simple case of parabolic dispersion, where we can 
do the calculation in two different ways: in the Landau 
level (LL) basis and in momentum (k) space. We then 
use the k-space approach to shed light on the crucial role 
of the dynamics of phase fluctuations. Finally we exam- 
ine quantum oscillations for arbitrary dispersion using 
semi-classical quantization of k-space results. 

Landau level analysis: Consider electrons with a 
dispersion et = h?h'^/2m* and chemical potential /i 
in a magnetic field H. In the LL basis Go{n,iuJi) = 
{ioji — Cra)^^, where w; = {21 + 1)ttT is the fermionic Mat- 
subara frequency and the spectrum C„ = (n+I/2)^c^M, 
with cyclotron frequency lOc = {eH/m*c) and LL index 
n. Using this Gq and the fluctuation propagator D^^, 
we obtain the self-energy (sec Methods) making the ap- 
proximation of retaining only diagonal terms in the LL 
sclf-encrgy matrix 'E{n,n';iiOi). Wc then calculate the 
DOS N{uj) = -{1/ttF) J2n ImG(n, itui ^ lo + iO+) with 
Lo measured from fi. All results shown here are for T = 0. 

We show in FigH^a) and (b) the DOS Nh{0) in a 
fleld H, normalized to the H = normal state No{0). 
Both curves in panel (a) exhibit quantum oscillations pe- 
riodic in {1/H) with the usual Onsager frequency given 
by the Fermi surface area. The index along x-axis di- 
rectly counts the number of filled LL n-p ^ fi/hwc- There 
is however a striking difference between the damping of 
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FIG. 2: Spectral functions: False color plots of j4(k, oj). a 
yl(k,aj = 0), for static phase fluctuations contributing to the 
self-energy E(k,0), showing intensity along Fermi arcs near 
the nodes. We use a parabolic dispersion with H — 0.003 (as 
in Fig. [1] we express H in units of /i = 1 and A ~ 0.3). b: 
Lj VS. k dispersion near the antinode (AN) along a cut per- 
pendicular to the FS, showing a Bogoliubov-like dispersion 
with a AN pseudogap. c: Gapless dispersion near the node 
(NN). We show in panel (a) both the AN and NN momentum 
cuts, d: A(k,iu = 0), for dynamic phase fluctuations (with 
r — v-p/t), shows the restoration of the full Fermi surface, e, 
f: Similar to plots (c,d) for the dynamical case. Note particu- 
larly the appearance of a zero-energy quasiparticle inside the 
pseudogap in the AN cut (panel e) due to quantum motional 
narrowing (see text). 



the static (black curve) and dynamic (red curve) results 
in Figllja). In the static case, the oscillations decay 
rapidly as exp (— tt/tWc) from the large intrinsic damping 
l/r ~ |ImE(nF,0)| 7^ arising from scattering of elec- 
trons from static phase fluctuations. (We note that the 
static c?-wave oscillations, though strongly damped, are 
still are much less so than the static s-wave results 13lll4| 
due to the averaging over the sign changes in the local or- 
der parameter.) In contrast, there is no intrinsic damping 
in the dynamic case, with ImI](7iF,0) = 0. The damp- 
ing in Fig[lja) arises from a small impurity broadening 
1/t ~ 7o put in by hand in Go, and inevitably present 
in real materials. Below, we will gain insight into why 
the intrinsic damping due to dynamic phase fluctuations 
vanishes. This has direct implication for quantum os- 
cillations in cupratcs, where no perceptible damping, in 
addition to that expected from impurity scattering, has 
been observed. 

In Fig. mb), we plot Nh{Q)/Nq{0) as a function of H 
(rather than l/H) for low fields with np of order hun- 
dred. Quantum oscillations are seen only for F ^ and 
completely suppressed for the static case. We also see 
a large, 77-dependcnt DOS suppression relative to the 
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FIG. 3: Fermi surface reconstruction a: Open and closed 
FS segments shown in the original Brillouin zone (BZ ) for 
density wave order scenario of ref. [l^ . We choose parameters 
t' /t = -0.4, 4>N = 0.143, V = O.llt and V' = 0.09t (enhanced 
l/'s to make the reconstruction visible) with ^ = — 1.05t. b: 
DOS quantum oscillations obtained from semiclassical quan- 
tization of the reconstructed FS in a vortex liquid state. Here 
t' /t = -0.4, (^iv = 0.143, V = 0.02f, V = O.Olt with 
/i = — 1.05t, so that the hole doping x ~ 0.12 and the closed 
orbit area ~ 2% of the BZ. The DOS also includes the con- 
tributions of open orbits. The _ff-dependent suppression of 
the DOS arises from phase fluctuations with A ~ 0.3t and 
r = 0.1i;f/^. Here 70 = O.Olt 



zero-field normal state, with, as we show below, a v H 
singularity at small H. Quantitatively, the suppression 
depends on F, becoming larger with decreasing F and 
most pronounced in the static case. The DOS suppres- 
sion also depends on the pairing strength A, which we 
take to be iJ-indcpendcnt, as is reasonable for low fields. 

Momentum space analysis: To gain insight into the 
LL results, we turn to a k-space analysis. The self-energy 
S(k, io;;) is calculated using the fluctuation propagator 
dll) with G'o(k, iw/) = {iuji — + m)"^! (see Methods 
for details). We first look at = /i^k^/2m* and then 
generalize to arbitrary dispersion later. 

It is helpful to look at the one-electron spectral func- 
tion A{'k,Lu) = -(l/7r)ImG(^)(k,a;). We show that dy- 
namical phase fluctuations restore a zero-energy quasi- 
particle [l^,!!^] at the antinodal kp. We may think of this 
as quantum motional narrowing, with the effect of pair- 
ing on the spectral function washed out on the longest 
time scales, as we now describe in detail. 

We plot in Fig. [5] (a,b,c) the spectral functions for 
static case with time-independent phase fluctuations, 
that should be contrasted with the corresponding results 
in panels (d,e,f) for the dynamic case with F 7^ 0. In 
the static case, we see in panels (a,b,c) that phase fluc- 
tuations broaden the node of the d-wave SC into "Fermi 
arc" , a region of gapless excitations where the spatial 
fluctuation- induced line width vp/^ overwhelms the gap 
Ak = A[cos(fc^a) — cos{kya)]/2. There is a pseudogap in 
the antinodal region where |Ak| > w_f/^ as seen in both 
panels . We note that the static case is not the (highly 
singular) limit F — at T = 0, but more closely re- 
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lated to the high temperature regime where T > F; see 
Methods. In the high-T regime, the phase fluctuations 
are classical |24| and we can ignore their time depen- 
dence. Thus the physics of the static results is relevant 
for high-T experiments like angle-resolved photoemission 
spectroscopy (ARPES). On the other hand, the quantum 
oscillation experiments are in the very different low-T 
limit, where the dynamics of phase fluctuations cannot 
be ignored. 

The results with F 7^ dynamics are qualitatively 
different from the static case. We see from Fig. [5] (d) 
that one recovers the full Fermi surface (FS), albeit with 
a highly anisotropic self-energy, as illustrated by the 
A{]<.,uj) dispersion plots in panels (e,f), along two repre- 
sentative momentum cuts perpendicular to the FS, one 
near the antinodc (AN) in the SC state and the other 
near the node (NN). We can understand the appearance 
of the zero-energy quasiparticle at the AN by looking 
at the self-energy. In contrast to the static case, which 
has a large antinodal |ImI](kF,w)| at low energies, dy- 
namical phase fluctuations lead to |ImE(kF, w)| ^ for 
<C F, the quantum motional narrowing mentioned 
above. The corresponding ImI](kF,a;) then leads to a 
quasiparticle pole at the chemical potential, even though 
the self-energy effects are are strongly k-dependent as 
seen from Fig. [2]|^e,f). 

The existence of sharp quasiparticles all around the 
full FS immediately leads to the quantum oscillations 
with the Onsager frequency. We define a renormalized 
dispersion Sk = Ck + S(k, 0) for low-energy quasiparti- 
cles, which has a non trivial 7J-dependence from the self- 
energy. We then use a semiclassical prescription [2^ to 
quantize the orbits (see Methods). The resulting DOS 
from this k-space analysis is shown in Fig. [Ijc), with 
a small impurity scattering 70 that damps the quantum 
oscillations. 

The most non-trivial aspect of this result is that 
the quantum oscillations ride on top of a large, field- 
dependent suppression of the DOS Nh{0), just as we saw 
in the LL analysis (Fig[ljb)). We can analyze this sup- 
pression by looking at the "average" DOS (without any 
semiclassical quantization), shown as the dashed curve 
in Fig. [Ijc). We can show analytically that the H- 
dcpcndent self-energy I](k, 0) leads to Nh{0) oc \/H as 
H ^ (see Methods) in the static case. This reproduces 
the celebrated Volovik result @ from quite a different 
route. The residual value of Nh{0) at = is due to 
the impurity scattering 70. 

Fermi surface reconstruction by a competing 
order: The analysis presented above shows that while 
phase fluctuations are able to reconcile quantum oscilla- 
tions with a large suppression of the DOS that goes like 
VH, they do not affect the oscillation frequency. Thus, 
to get a complete description of the underdoped cuprate 
experiments, we need to incorporate both phase fluctu- 
ations and a competing order. We can incorporate any 



one of the proposed broken symmetries 17H2ll] within 
our k-space formulation. Only experiments will decide 
which competing order is most relevant for a particu- 
lar material. Once SC long-range order is destroyed by 
the field, it is natural that the ground state of a lightly 
doped Mott insulator exhibits a new density wave insta- 
bility that reconstructs the FS. However, we are firmly 
of the opinion that the large (~ 50meV) AN pseudogap 
cannot arise from a small (or subtle) symmetry breaking, 
and it is not reasonable to use a large symmetry breaking 
potential to reconstruct the FS. 

To understand the interplay of FS reconstruction and 
phase fluctuations, we analyze, as an illustrative ex- 
ample, the density- wave order proposed in ref. [l^ 
with Q = (f ,0) consistent with recent high- field NMR 
data [ig Following ref. [l^, we start with Ck = — 2t[(l -|- 
(/)n) cos kxa + {l — (/)N) cos kya] — At' cos k^a cos kyU, where 
t {t') is the nearest (next-nearest) neighbour hopping and 
(j)N the nematicity. Phase fluctuations renormalize this 
dispersion to ek = Ck + S(k, 0), with the self-energy dis- 
cussed above. To find the FS reconstruction, we diago- 
nalize the Hamiltonian 



Hk 



V 

v 



V 
V 

v 



V V 

V v 

ek+2Q V 

V ek+3Q / 



where V and V are the density wave potentials. 

We semiclassically quantize the resulting energy dis- 
persion. The results are shown in FiglS) The chemical 
potential is fixed for x w 0.12. As shown in Figl3]a, there 
are small electron-like FS pockets and open FS sheets. 
The pockets have an area around 2% of the original BZ 
and there is only one pocket in the new BZ. Although the 
open FS sheets do not contribute to oscillations of DOS, 
they give rise to large contribution to DOS in the absence 
of pairing. In a vortex liquid these DOS from the open 
orbits are largely suppressed as shown in Figl3]b. Hence 
in this case the frequency as well as the suppression of 
DOS are both in good agreement with experiments [iHl]. 

In conclusion, we have presented a calculation of the 
electronic excitations in a field-induced vortex liquid 
state with short-ranged d-wave pairing fiuctuations. This 
leads to a non-trivial self energy that is responsible for a 
large suppression of the DOS, with a constant plus ^/H 
variation at low fields. The dynamics of phase fluctua- 
tions restores the zero-energy quasiparticles via quantum 
motional narrowing at low temperatures, thus leading to 
the full FS, absent in the zero-fleld, high-T normal state. 
A competing order parameter that breaks translational 
symmetry then reconstructs the FS to give the observed 
low frequency of the quantum oscillations. 
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DE-SC0005035. 



5 



METHODS 

Dynamics of Fluctuations: We characterized the vor- 
tex hquid state by the ansatz ^ that describes short 
ranged d-wave phase fluctuations. The gaussian spatial 
decay is motivated by Ginzburg-Landau theory. The dy- 
namics of the fluctuations has the form 4- F) in 
Matsubara space. The static approximation, discussed 
in the text and used previously for s-wave SC's [lil. [l^. 
corresponds to {l/T)Sui,.o. The static case is thus related 
to the high temperature regime T ^ T, and not the 
(highly singular) F — > limit of the dynamics case. At 
high T, the classical, thermal fluctuations that dominate 
are time- independent. 

h-space analysis: We find the self-energy in the static 
case 



S(k,a;) 



2 Vi, 



exp 



2vl 



(2) 



where Fj:){x) = —i^yTTeri{ix)/2 is the Dawson function. 
Here Ak = A(cosfc2;a — cos/cya)/2, (,k = ek — fj. and Vk = 
|9^k/9k|. Derivations are sketched in the Supplementary 
Information. 

In the dynamic case, the imaginary part of E is 



Iml](k, uj) = 
tami 



Alt 



2T 



2'KVk 

n 



f]2 + r2 



exp 



coth 
(a. 



2f) 



H 



(3) 



We take the T = limit of this result and obtain the real 
part by the Kramers-Kronig transform of ([3|) . 

\/H-hehavior: At sufficiently low fields {i — > oo) we 
may use the static approximation. As H — > 0, the 
static self-energy ([2]) reduces to the d-wave BCS-Gor'kov 
form Y.(k,iuji) ~ A^/(za;; + ^k)- For small H we ex- 
pand E(k,w = 0) in p = k — kjv around one of the 
four nodes at k^r, with ~ P-va and p± ~ p.v-p. 
We get ReE(k,0) w {vl/vF)plp±i^ and ImE(k,0) w 

and v-p are evalu- 



-y^{v\/vF)ple, where va 



I 9Ak 
9k 



ated at k^r. Using this field dependence of E(k, 0) in the 
spectral function leads to Nh{u} = 0) oc \/H as H 0, 
thus recovering Volovik's result Q. 

We have also generalized the above asymptotic analysis 
to include impurity scattering 70. To make the algebra 
tractable, we use a Lorentzian form for ImE, rather than 
the Gaussian in ([5]). We then find a result of the form 
Nh{0) = Ajo + B\/H as — > 0, in excellent agreement 
with the numerical results shown in the text. 

Semiclassical quantization: We use the semiclassical 
prescription for quantizing electron orbits with the renor- 
malized dispersion ek = ek + S(k, 0). For a given field H, 
we generate a set of energy levels {?„} as the solutions 



of A(en) — {2Trn/£'^), where A{T) is the k-space area en- 
closed by a closed orbit at energy and n is a positive 
integer. We then use the ?„'s to compute the DOS. 

LL Analysis: The d-wave self-energy T,{n,n';iLUi) is 
a matrix in the LL basis, in contrast with the s-wave 
case [1^. For large LL index n, the matrix elements de- 
cay like exp[— (n — n')'^/An] away from the diagonal. In 
the spirit of our calculation, which does not treat fluctu- 
ations self-consistently, we neglect off-diagonal terms and 
retain only the diagonal ones, denoted by E (71,10;;). For 
the static case, the self-energy at large n is given by 



E(n,c.) = — 



X 1 + Jo (2Vn + ni^) - 2Jo (^^/2{n + m)^'^ 



(4) 



where Inm = (n + ni)!/n!rti!2"+"i+^ and Joix) the ze- 
roth Bessel function. 
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SUPPLEMENTARY INFORMATION 
Electronic Green's function in a d-wave vortex liquid 

Here we give some technical details of the calculation of electronic self-energy and Green's function in a vortex 
liquid. The d-wave pairing is described by the field 4'^(r -I- a/i/2,T) defined on the link (r, r -I- ajj,). Here r is the 
imaginary time with < r < /3 = l/T, the inverse temperature (We set k-Q = Ti ^ 1). The single particle Green's 
function G and anomalous Green's function F satisfy the Gor'kov equations: 



G(r,T;r',r') = Go(r,T;r',T') - ^ j ^drirfrs j dvidv2 



Go(r,r;ri,Ti)*^(ri + ^,ti) 



Go(r2 -Fai>,T2;ri + a/i, n)** (r2 + — , r2)G(r2, T2; r', r') 
F(r,r';T,T') = j driGo(r, r; n, ri)4'^(ri + ^, Ti)G(r', r'; n + a^, r') 



(Sla) 
(Sib) 



The indices /x, v run over bond-directions ±i,±y. The factors 1/16 and 1/4 are related to the normalization of 4*. 

The self-energy shown in the inset of Fig. 1(a) of the paper is obtained as follows. Generalizing the s-wavc approach 
of ref. [H, we first average the above equations over the configurations {^'^(r,T)} and then make the decoupling 
approximation 

(*^(ri + ^,r0*:(r2 + y , r2)G(r2, r2; r', r')) « (*^(ri + ^,ri)*:(r2 + ^, r2)) (G(r2, r2; r', r')) (S2a) 



(*^(ri + ^,Ti)G(r',T';ri +aA,ri)) « (*^(ri + ^, n)) (G(r', r'; n + ^A, n)) 



(S2b) 



By definition (^'^(r, r)) = in a vortex liquid, which implies that (F(r, r; r', r')) = 0. 

For notational simplicity, we denote (G(r, r; r', r')) by G(r, r; r', r'). Eq. (jSla^ can then be written in the form 



G(r, r'; iuji) = Go(r, r'; iuJi) - X! X! ^ 



Dau r2 - ri 



16 

au ~ aji 





, 2e 


1 


exp[— i — / 







A • dl]Go(r,ri;iwz) 



i^lm Go(r2 + aO, ri + afi; -iuji + iri„i)G(r2, r'; i^z) 



(S3) 
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Here tui = {21 + l)7rT and fli = 21ttT are Fermi and Bose Matsubara frequencies, respectively and J^^ denotes J dri. 
Using the the separable form of the gauge invariant fluctuation propagator (Eq. (1) of the paper), we find that 

D^,{r, i^rn) = s^,2?(r)^(^O„0 = s^^A^ cxp (-rV2£') | (S4) 

Note that J^(irim) is the Matsubara representation of the dissipative form ImJ^^^^n) = + F^). 

The electromagnetic phase factor appearing in Eq. (|S3p results from the definition of gauge-invariant propagator 
defined in main text. The integral A.dl is taken along a straight line connecting and r. The reasoning for 
this parallels that given in ref. [l[ for the s-wave case. One additional complication here is the phase factor = 
{a/2i)^{vx + ^x)(yy — fly) arising from the definition of the d-wave pair field on the bonds of the lattice. (That this 
is a lattice related effect is evident from the fact that Q^^, is negligible in the low field limit, a <C i-) 

We choose to work in the Landau gauge A = Hxy to represent the magnetic field Hz. For a parabolic dispersion 
the Green's functions in the Landau level (LL) basis is G{nq,n'q';iuji) = / o?rc?r'(/)*g(r)^„/g/ (r')G(r, r'; iw;). Here 
0„g(r) is the LL wave function with n the LL index and q going over the degenerate states in each LL. We can 
now rewrite Ea. (jS3[) in the LL basis so that it has the form of Dyson's equation G = Gq + GqSG. The bare 
Goinq,n'q';iLJi) = Go{n,iuji)Snn'Sqq' with Go{n,iuji) = {iuji - ^„)"^ The sclf-encrgy is given by 

S(ng,n'g';ia;,) = - ( — ) ^ ^ ^ s^i.V(n(j, nV, ni)^ _ " . \ (S5) 



Ii_i,y{nq,n'q\ni) 



E 

91 



exp[— I — 



A • d\\D r2 - ri 



ri + - 



av 



, (r2 + ai>)0;^g^ (ri + afl)(j>n'q, (ra) 



The Matsubara sum can be done in the dynamic case (F 0), using the form of the J-{iQrn) in Ea. (jS4[) . In the static 
approximation £)^i,(r,r) is independent of r so that J^{iilrn) = Sm,o/T. As a result, the Matsubara sum in Eq. ([55)) 
is trivial. We discuss below how the static limit is recovered by taking a suitable high-temperature or classical limit 
of the dynamic case. 

We can do the J2qi r2 integrals and express /^y (ng, n'q\ ni) in terms of special functions, with the various phase 
factors "canceling out" . We omit the details of this lengthy algebra here. The self-energy evaluated from Eq. (|S5[) 
turns out to be diagonal in g-space. It is also independent of q, which can be traced to the fact that the vortex-liquid 
state does not break translational invariance. We denote it by 'E{n,n';iuji) in the main text. Unlike the s-wave 
case [ij, E is a matrix in the LL index due to the d-wave nature of the order parameter. We can show that the matrix 



elements n'; iuji) decay rapidly away from the diagonal n = n' like ^ exp [—{n 



,l\2 



/4n]. In the spirit of our 



non-self-consistent calculation [Fig. la (inset)], given that Go is diagonal in the LL index, we only retain the dominant 
diagonal terms in S(n, n; iuii) = S(?t., iuJi). The k-space analysis described below serves to validate this approximation. 
In the dynamic case, the imaginary part of the retarded self-energy (iw; — > w -I- iO"^) is 



ImE(n, a;) 



A 
2" 

2^ 



E 

coth 



2e" 



a 

2P 



2T 



tanh 



2T 



ImGo(ni, r2 - a;)ImJ"(-") (O) I 



(S6) 



Here Inni = {n + ni)!/n!ni!2"'+"i+^ and Ln{x) is a Laguerre polynomial. The real part of the self-energy can be 
obtained using Kramers-Kronig relation from Iml](n, w). For large values of the LL index n we can derive a useful 
approximation 

2 

r / , n\ / n\l I I / 



ImS(7i,w) - - ( y 



l + Jo {2V2^^^ -2 Jo (2\^^ 



4:\/mrhiiJr 



coth — 
2T 



— tanh 



2T 



n 



(S7) 



where Jq (x) denotes a Bessel function. We have benchmarked this form by comparing it with results obtained directly 
from Ea. ((S6| and find that the approximation is accurate for up — {fi/fnoc) > 10. We have used Eg. ijST)) for computing 
the results reported in the main paper. We have also found expressions analogous to Eg. dSSI) and Ea. (|S7[) for the case 
of static phase fluctuations separately, and summarized in the Methods section. 
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k-space calculation 

We now discuss the k-space calculation of self-energy I](k, w) used in the main text together with a semiclassical 
approach to quantum oscillations for /i > A 3> huc- As a result of the translational invariance of the vortex-liquid 
state, one can make a gauge transformation and represent the (non-gauge-invariant) Green's function G(r,r']iuji) as 



G(r, r'; iw;) = G(r' — r, ioj;) cxp[— i- / A.dl], 



(S8) 



where G depends only on the separation r' — r. We rewrite Eq.(|S3P in terms of G{r,iuji) or its Fourier transform 
G(k, loj;) (where wc omit the tilde to simplify notation) and obtain the corresponding sclf-cncrgy 



E(k, iioi) = -T 



■ ^X'(q)J'(ifi„i)Go(k - q, ~iuji + iO„i) 





, qxa\ 




cos ^ 


'^xCl — 

2 / 


— cos 



2 



(S9) 



Here I?(q) = cxp(— q^^^/2) is the Fourier transform of I'(r) in Eci.([T]). 

In the semiclassical regime where jj, ^ hojc or the cyclotron radius £ the inter- vortex separation, we neglect the 

effects of quantization of electronic orbits while evaluating S(k, iw;) and approximate Go(k, iuji) by (iw/ — ^k)~^. We 
can further simplify this result by noting that in the field regime of interest £ 3> ^-nd I'(q) is sharply peaked around 
q = with a width ~ l/£. This enables us to make the expansions cos {k^a — ^|^) ~ cos(fc3,a) and ^k-q — ~ v^.q, 
where Vk = (9^k/i9k). Proceeding in a manner similar to the LL discussion above, we obtain for dynamical fluctuation 
case 



ImS(k,cj) = -- ^ 



2V27rvk 



dn 



coth — — tanh 

2T \ 2T 



fi2 + r2 



(SIO) 



Here Ak = (A/2)(cos fc^ja — coskya) and I'k = |vk|- The real part of S](k,a;) can be obtained by Kramers-Kronig 
relation. The corresponding expression for I](k, cj) in the static case is shown in the Methods. 



The high-temperature (classical) limit and the static case 

As noted in the paper, one cannot just take the (singular) F — > limit of the dynamic phase fluctuations at low 
temperatures to describe static phase fluctuations. Here we show the relation between the high-temperature limit of 
the self-energy of Ea. (|S10p for the dynamic case and static case self-energy of Eq. (2) of the Methods section. In the 
high temperature limit T ^ T, we see that T ^ for the frequencies contributing to the integral of Ea. (|S10p . Using 
coth(rj/2T) ~ 2T/n and tanh(rj/2T) ~ n/2T we get 

TAli f 1 



TTUk 



ImS(k,w)c.^=^ / dl^^^^^e (Sll) 



To see that this is qualitatively similar to the self-energy for the static case wc approximate the Lorentzian 1/(^1^ + r'^) 
above by a Gaussian ^7r/2exp(— 172/2F2)/F2. We thus obtain 

T^y-(k^,„) ^ ( ^ ) J- —-^ ^ 2(„^r^.^) ^S12) 

^ ' \ r jy 2 + r^£2 ^ ' 

This has exactly the same form as the static case (Eq. (2) of Methods) with some redeflnitions: (TA^/F) going to A^ 
and V^'k + r^^'^ to Uk. For a closely related discussion, see rcf. 



DOS in the low field regime: \^ plus impurity scattering 

To investigate the asymptotic behavior of the DOS, Nh{0) as a function of flcld and impurity scattering strength, 
wc focus on the self-energy for the static case. To obtain analytical results, wc approximate the Gaussian in Eq. (4) 
of Methods by a Lorentzian, so that we get 
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S(k,a;) 



a; + + j(7o + I'k/-^) 



«7o 



(S13) 



where 70 <C A < // comes from weak impurity scattering. This corresponds to a static pair correlator X'('r) 
exp(— r/£) (that decays exponentiahy rather than as a Gaussian) and gives rise to qualitatively similar low-energy 
electronic spectra and DOS. We can also use a more general form of (|S13p with 7(cj), but only its w = limit 70 is 
relevant below. 

We calculate Nu{0) as 



iV(0) = - 



1 

2^3 



(fk 



ho + VI./ em 



+ 7o' + 70^'k/^) - e^k/^ 



(S14) 



We next expand Ak ~ vaP\\ and ~ vpp±, with p = k — k^r, in the same way near the nodes kjv's as described 
in Methods (where we did not include impurity scattering). In the low- field regime H ^ 0, we obtain Nh{0) by 
expanding in powers of (l/i) to get 



iVff(O) 



70 



dp±dp» 



VF 



I 2 2 



(S15) 



where the momentum cut-off po — 7r/a. This implies that Nh{0) 
appropriate constants that depend logarithmically on jq. 



A'jo + B\/H as if — > 0, where A and B are 
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